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£3 1 Abstract 
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For the color-suppressed B° — ► D°ir° decay, nonfactorizable contributions are expected to be 
leading and naive factorization description breaks down. We study l/m& power-suppressed non- 
factorizable effect in B° — ► D°-k°, which is due to soft exchange between the emitted heavy- light 
quark pair and Bn system, in the framework of QCD light-cone sum rules (LCSR). The resulting 

<D : 

correction to the decay amplitude is found to be numerically comparable with the corresponding 
factorizable piece, estimated at about (50 — 110)% of the latter. The relevant parameter 02 receives 
a positive number contribution, due to the factorizable correction and the power-suppressed soft 



effect. Our findings would be crucial to phenomenological understanding of the B° — ► D°ir° 
decay. 
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I. INTRODUCTION 

"Naive" factorization []], or "generalized" factorization [3] developed subsequently, had 
been viewed as a simple but predictive model for two-body hadronic decay of heavy mesons 
prior to presentation of QCD factorization Q]. At present, it is known to us that for a 
large class but not all of two-body nonleptionic B decays, QCD factorization can furnish 
a rigorous theoretical basis for the naive factorization assumption of the hadronic matrix 
elements. Some of examples, for which the naive factorization holds up to power corrections 
in AgcD/n^b and a s , are the charmless decays B — > tttt, ttK and the class- 1 charmed decays 
B° — > D^ + n~ (relevant to the parameter a\). In the heavy quark limit m& — > oo, but 
to all orders of perturbation theory, this type of processes can be systematically computed 
in terms of convolutions of hard-scattering kernels with the leading twist light-cone distri- 
bution amplitudes of the corresponding emitted light mesons. Such a treatment is based 
on color-transparency argument that for the limit in question the momentum carried 
by an light meson which is directly emitted off the relevant weak vertex is so large that it 
has not sufficient time to exchange soft gluons with the system including the decaying B 
meson and the produced meson picking up a spectator quark. Typical examples for which 
QCD factorization does not apply are the class-2 charmed decays B° — > D^°ti (usually 
called color-suppressed decays for the reason that the relevant phenomenological parameters 
a2{Dn) is of O(l/N ) in the large N c accounting |5[). The reason is that unlike the aforemen- 
tioned case where a light meson is emitted, the heavy meson, as an emitted particle 
which is neither small ( ~ 1/Aq C d) nor fast, is produced in the color-suppressed decays and 
so can not be decoupled from the Bw system. This indicates clearly that the factorization 
contributions to B° — > D^tt don't provide a leading result, for all that they are present, 
and nonfactorizable soft contributions, for example, the charge-exchanging rescattering pro- 
cesses p from the dominate class-1 channel, dominate in such decays. Hence, at the present 
a theoretical understanding is not accessible for the color-suppressed decays B° — > D^°ir 
and also for the class-3 decays, say, B~ — > D^°tt~ to a certain extend, although in the 
latter case the class-2 amplitudes are predicted to be power-suppressed with respective to 
the corresponding class-1 ones in QCD factorization. It is quite a challenge to give a consis- 
tent theoretical explanation of the data on B — ■> D^ir, after new experimental observations 
B{B° -> D°7T°) = (2.74^;| ± 0.55) • 10" 4 7] and B(B° -> D ^ ) = (3.1 ± 0.4 ± 0.5) • 10" 4 
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are announced respectively by the CLEO and Belle Collaborations. While there exist 
the attempts to quantitatively understand the color-suppressed decays B° — > D^ %° in 



perturbative QCD(pQCD) toll and soft-collinear effective theory (SCET) [lOj], the model- 
independent discussions jlll llj. can help to get interesting suggestions about the mag- 
nitudes and relative phase of the parameters ai )2 . A couple of important findings can be 
summarized as follows: (1) A sizable relative strong interaction phase is expected between 
class- 1 and class-2 B — ► D*n decay amplitudes (2) The parameter |a 2 | is extracted to 
be |a 2 (L>7r)| ~ 0.35-0.60 and \a 2 {D*ir)\ ~ 0.25-0.50 from the data [ljj (3) Several types 
of possible power corrections to a± parameter have been estimated and found to be small, a 
near-universal value |ai| ~ 1.1 observed experimentally is now put on a firm footing jll]. 

Now we are not able to give a reliable theoretical interpretation for the first two obser- 
vations, because of the unknown leading nonperturbative effects involved in the parameter 
a 2 (-D7r). However, the l/m b power-suppressed nonfactorizable contributions to a 2 (L>7r), 
which come from soft exchange between the emitted heavy-light quark pair and the Bit 
system, would be expected to be much more important than in the case of B — > ttti, 7rK and 
thus a reliable estimate of such effect is crucial. 

Earlier discussion on the power-suppressed contribution to B° — > D°ir in the QCD light- 
cone sum rule(LCSR) approach [l^ | can be found in Ref. Q]. Also, in the framework of a 



generalized QCD LCSR 
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the similar effects have been estimated for some of the other 



important B decays 12, Q|. In this paper, we intend to apply the generalized QCD 
LCSR approach to estimate the soft effect in the color-suppressed B° — > D°n° decay and 
then compare the result yielded with the naive factorization contribution. 

This paper is organized as follows: the following Section contains a detailed derivation 
of LCSR's for the power-suppressed soft contribution to the B° — > D°ir decay amplitude 
and the numerical results. The last Section is devoted to discussion and conclusion. 

II. LSCR'S FOR SOFT NONFACTORIZABLE EFFECT 



The relevant effective weak Hamiltonian for the B — > D ir decay is written as 

Gp 

71 



n w = % v^Ac^o^) + c 2 (/i)o 2 (/i)l , (i) 
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where are the Wilson coefficients, the CKM matrix elements and the four-quark 
operators given by 

o x = (cr^)(Jr M 6), o 2 = (dr"u)(cTpb), (2) 

with = 7^(1 — 75). Further, the use of Fierz transformation can make ([!} rewritten as 



G 



v/2 

where 



{CM + iCaO*))©!^) + 2C7 a (A*)6i0i)] , (3) 



Ox = (c^r,u)(d^m) (4) 



with A a being the color 577(3) matrices. 

Among the non-leading part of the decay amplitude (7} 7r°|77w|-B°) is the factorizable 
and power-suppressed soft contribution. We can parameterize it as 

A NL (B° — 7JV) = A F (B° — > £>V) + A(5° — 7JV) (5) 
= -^^K>l/D^(m^ L . 

where F B7r is the B — > n form factor, /d the D meson decay constant, Af{B° — > 7J°7r°) 
and ^4.5(5° — ► D°n°) express the factorizable and power-suppressed soft contributions 
respectively: 

A F (B° — 7JV) = ^^JfK)^ + ^), (6) 

^(75° — > 7JV) = ^V cfe \C,(2C 2 )(7JV|d 1 |75 ) 5 , (7) 
and the parameter is defined as 

For a quantitative estimate of the nonfactorizable matrix element (7J°7r°|Oi|77 ), we make 
use of the generalized QCD LCSR method developed in [la] . We start with the correlation 
function: 

F a (p,q,k) =i 2 [ <fxe~^ x [ d i ye i ^(7r ( q )\T{j^(y)d 1 (0)ji(x)}\0) , (9) 



where = W7a75c and jf = m^bi^^d are currents interpolating the D° and 77° meson 
fields respectively. The correlator is a function of three independent momenta chosen to be 
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q, p — k and k. An important point of this method is to introduce a fictitious unphysical 
momentum k. Consequently, in the correlator the quark states before and after the 6-quark 
decay will have different four-momentum, and thus one avoids a continuum of light parasitic 
contributions in the B channel. Of course, the unphysical quantity must disappear from the 
B° — > D°tt° ground state contribution in the dispersion integral. This can be guaranteed, as 
will be seen, by picking out a reasonable kinematical region for which the LCSR calculation 
is effective. 

The kinematical decomposition of the correlation function © can proceed in the following 
form: 

F a = {p - k) a F^ + q a F® + k a FM + e a(3Xp qPp x k?F^ , (10) 

Here, F l, s are scalar functions of 6 independent Lorentz invariants, which are chosen to be: 
P 2 = (p + q — k) 2 ,p 2 ,q 2 , (p + q) 2 , k 2 and (p — k) 2 . Using the Operator Product Expansion 
(OPE) near the light-cone x 2 ~ y 2 ~ (x — y) 2 ~ 0, the correlator Q is calculable. For the 
calculation to go effectively, the momenta squared P 2 , {p + q) 2 and (p — k) 2 have to be taken 
spacelike and large in order to stay far away from the hadronic thresholds in both the B and 
D channels. Furthermore, a simple and possible choice, for the external momentum squared 
k 2 and kinematical invariant p 2 , is to let k 2 = and p 2 = m 2 D , tub being the mass of D°. 
The pion is taken on shell and we set q 2 = 0. Altogether, the kinematical region used for 
our LCSR calculation is 

q 2 = k 2 = 0, p 2 = m 2 D , \(p - k)\ 2 > A QCD , \{p + q)\ 2 > A QCD , \P\ 2 > A QC d • (11) 

In this region the light-cone OPE is applicable to the correlator (0) and the result can 
be expressed in the form of hard scattering amplitudes convoluted with the pion light-cone 
distribution amplitudes. We note that in <fT0|t the relevant invariant amplitude to our desire 
is only F^ p ~ k \ The QCD result for F^ p ~ kS) is, in a general form of dispersion relation, 
expressed as 

j?<p-k),, M 2 / , ^2 o2n 1 lm s F^cp(s,(p + q) 2 ,P 2 ,p 2 ) f . 

F$ CD \(P -k),(p + q),P) = -J mi ds s _ {p _ k y • (12) 

On the other hand, we can obtain a corresponding dispersion relation in the hadronic 
level. By inserting in the right hand side of a complete set of hadronic states with D° 
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quantum numbers, we get: 
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where p^(s, (p + <?) 2 , P 2 ,p 2 ) and are respectively the spectral function and the threshold 
mass squared of the excited and continuum states in D channel, H((p + q) 2 ,P 2 ,p 2 ) is a 
two-point correlation function with the following definition: 

U((p + q) 2 ,P 2 ,p 2 )=i J d 4 xe-^ + ^(D\p-k)n%q)\T{d l (0)ji B \x)}\0). (14) 

By assuming quark-hadron duality we replace with the effective threshold of the 
perturbative continuum Sq, and substitute the hadronic spectral density p% in |T3|l with the 
corresponding QCD one, i.e : 

p°(s, (p + g) 2 , P 2 ,p 2 )6( S - O = hm s F%c$(s, (p + g) 2 , P 2 , p 2 )9(s - tf) . (15) 

Matching the hadronic relation (jlHj) onto the QCD result ()12|1 yields the expression 

^n(( P + g) 2 ,p 2 ,p 2 ) _ i ^^ to^jp+j^ (16) 



m D ~~ (P ~~ ^) 2 71 ^ s ~ (P ~~ ^) 2 

which then becomes 

H((p + g) 2 , P 2 , P 2 ) = -4" f J ^ e ( ^- s)/M2 Im s pfe fc) ( S) (p + q) 2 , P 2 ,p 2 ) , (17) 

after the Borel transformation in variable (p — k) 2 . 

Next, for the above expression which is only valid at large spacellike P 2 , we have to 
perform an analytic continuation to large timelike P 2 , keeping the variable {p + q) 2 fixed. A 
natural continuation point is P 2 = m 2 B . The analytic continuation of 1)171 ) yields the result, 

U((p + q) 2 ,m 2 B ,p 2 ) = i J d'x e-^ + ^{D\p-k)Tx\q)\T{dMA B \m^) 

= f S " ds e^- s y M ^ s F^\s,(p + q ) 2 , m l,p 2 ). (18) 

Then we employ the analytical property of the amplitude U((p + q) 2 , m 2 3 ,p 2 ) in the spacelike 
variable {p + q) 2 - Inserting in the right hand side of (JTSj) a complete set of hadronic states 
with the B° meson quantum numbers, we have the following dispersion relation: 

m B - [P + Q) Js h s' - (p + q) 2 



where the B meson decay constant is defined as 

(B |6i T5 d|0) = m B f B . (20) 

At this point, we would like to emphasize that the unphysical momentum k disappears from 
the ground state contribution due to the simultaneous conditions P 2 = (p + q — k) 2 = m 2 B 
and (p + q) 2 = m 2 B , so that the physical B° — > D tt° matrix element of operator 0\ is 
recovered. 

Then ({IB}) can be changed to a form of double dispersion relation as 

u(( P + q) 2 , m 2 B , P 2 ) = -4r f E dse^i M2 

x s '-{p + q) 2 l ™s> l ™A P CD\s,s>,m 2 B ,p 2 ). (21) 

The upper limit R of integration in s' rests generally on s, m 2 B and p 2 . At present, we make 
use of quark-hadron duality once more and approximate the integral in by the s' > s B 
part of the dispersion integral ()21|). where s B is the effective threshold in the B channel. 
After the Borel transformation with respect to the variable (p + q) 2 is made, LCSR for the 
B° — > D°7T° matrix element of the operator 0\ is of the following form: 

(D°(p-k)-K°(q) B°(p + q)}= ~] 2 f'° ds e^-V"* 

K 2 fDfBm B Jml 
x /' /?(S ' m| ' p2 ' SoS) ^„'.(- 2 R -')/M' 2 T^ .T^ p(P-fc) 



/ ' S ° W eK-^V, Im s F«( S , m|,/) , (22) 

o 

where R is the upper limit of the integration in s' after the use of duality ansatz. 

In order to obtain a LCSR estimate of the desired soft nonfactorizable contribution 
(D 7T \Oi\B°)s, which is due to soft gluon emission off the emitted heavy-light quark pair 
and subsequent absorption into the Bit system, we have to know the explicit expression of 
Fqcd ■ m the derivation of Fqcd we employ the light-cone expansion form with higher- 
twist terms included for a massive quark propagator, which, in the fixed-point gauge and 
only considering correction of operators with one gluon field, reads j^j], 

S ij (x 1 ,x 2 \m) = -i(0\T{q i {x 1 )q j (x 2 )}\0) 



X 



2 (k 2 — m 2 ) 2 k 2 — m? 



where G pu is the gluon-field strength, and g s the strong coupling constant. This means that 
only the higher-twist components of light-cone distribution amplitudes for the relevant pion, 
which are corresponding to quark-antiquark-gluon nonlocal operators, would be involved in 
the final LCSR result. 

After some lengthy calculation we obtain the twist-3 contributions: 



4v / 2tt 2 •'O J % m 2 — (p + q(l — a{)) 2 Jo (1 — x)m 2 c — Q 2 x(l — x) 
x q ■ (p — k) (2 — v)q ■ k + 2(1 — v)q ■ (p — k) 

mbh - r x dvlvcu n t**^ f 1 dx x{1 ~ x){2x ~ 1] 



A\f2-K 2 Jo J 1 m 2 . - (p + q(l - ax)) 2 Jo m 2 c - Q 2 x(l - x) 
x q ■ (p - k) [(2 - 3v)q ■ k + 2(1 - v)q ■ (p - k)] . (24) 

where (p^ is a twist-3 distribution amplitude of pion. The definition of the twist-3 distribu- 
tion amplitude as well as of the twist-4 ones (fi±, <p\\, (fi± and (f>\\, which will be encountered 
in the calculation, is given below through the relevant matrix elements: 

- \^(Q\d{Q)a^ b G a p{vy)d(x)\ix°{q)) = if 3lT [{q a q^gp v - qpq^av) 



-V2(0\d(0)i ltl G a p{vy)d(x)\7r°(q)) = q ^ QpXa f v [ Va^{a h ^-n^+v^) 

+(9iaHs-9ip<la)Uj Va^xia,, fi)e-^ + y^\ (26) 

-V2(0\d(0) lfMl5 G a p(vy)d(x)\n°(q)) = g /^^" /,/ Da^nh./ije-^ 1 ^' 

+{9taQfi - Sfo<l«)f* J Va t M^ ^)e~^ + y^ , (27) 
where / 37r is a nonperturbative quantity defined by the matrix element {^ua^v^G^d^), 
G a p = 7}e aj3pr7 G pcr , G p(J = g s X a /2G pa , Von = daida 2 da 3 5(l - a x - a 2 - a 3 ), and = 
g a p — (x a q/3 + X/3q a )/qx. The asymptotic forms of all these distribution amplitudes are given 



as 

, 2 



hn{oci,fi) = 360aia 2 a 3 , (28) 

(j>±(oa, fi) = 105 2 (/i)(«i - a 2 )a\ , (29) 

0||(ai,/i) = 1205 2 (n)e(fx)(ai — a 2 )aia 2 a3 , (30) 

4>j_(ai,ii) = 105 2 (/i)a 2 (l - a 3 ) , (31) 

0||(«i,/i) = -405 2 (/i)aia 2 a;3 . (32) 
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with <5 2 (p) and e(p) being two nonperturbative parameters. 

By changing the order of the integral variables, (|24|) is converted into the following form: 



7 {p-k) 
tw3 



™> b f: 



3?r 



ds 



1Q\/2tt' z Jml 
do 



dy 



du 



s — (p — k) 2 J^f- Jx(s,y,P 2 ) ml — (p + qu) 2 



X 



Jx(s,y,P 2 ) V 2 



c(s,y,P 2 ) 
+ 



3?r 



— it, u — t>, t>) s + ( (p + q) — p ) (2v — x(s, y } P )) 

y v 
i 



ds 



16v/2vr 2 Jml s-(p-k) 



2 m 



dy{2y-l) 



du 



x(s,y,P*) m\- (p + gtt)* 



X 



fit' 



rm 



z(s,y,P 2 ) IT 



— 03.(1 -u,u-v,v) -^-s+((p + q) 2 -p 2 )(2v-3x(s } y } P 2 )) , (33) 



where x = (s - ^)/(s - P 2 ). 

The derivation of twist-4 contributions turn out to be even more tedious. The result is 
given as follows 



ds 



, dy 



du 



x 



8\/27r 2 Jml s - (p - A;) 2 Jx(s, y ,P 2 ) ml - (p + gw) 2 

2\ 



— 0J_(1 — M, M — t>, f ) 

x(s,y,P 2 ) f 

"^6 A 



8V/27T 2 ^ s - (p - A;) 2 Vsi 



3 x{s,y,P 

v 

dy(2y-l) 



du 



(s,y,pz) m\- (p + g«) 2 



dt> 



x(s,y,P 2 ) V 

00 (is 



0_l(1 — n, m — f , f ) 



SV^TT 2 An? S - (p - k) 2 

dv ~ 

X / — $i(l 



2 



3 x(s,y,P 2 

V 

du 



x(s,y,P 2 ) V z 

00 ds 



x(s, y ,P 2 ) [(ml - (p + gw) 2 ]" 

2 



y 



■s - — + ((p + g) 2 - p 2 )(-w + x(s, y, P 2 )) 



x 



8V2~W m 2 s-(p-A;) 2 

di> 

c(s,y,P 2 ) V 2 

r°° ds 

„2 



dy{2y-l) 



du 



&\/2lX* Jiiif 



$i(l -u,v) 
1 



x(s, y ,P 2 ) [(ml — (p + gw) 2 ]^ 

2 



dy 



2(s ^) + (P 2 -s)t; 

1/ 

du 



s - (p - A;) 2 Jx(s, y ,P 2 ) [ml - (p + gw) 2 ]' 



x- 



$9.(V) 



s * + ((p + g) 2 - p 2 ) (-« + x(s, y, P 2 ) 



y 

ds 



8V2ir 2 Jm 2 [s - (p - k) 2 \ 



1 p2 



p2 m c 

y 



du 



x- 



$ 2 M 



x(s,y,P 2 
-P 2 



(2g-(p-A;)) 2 1 



(s,j/,p 2 ) m 2 - (p + gw) 2 
x(s, ?/, P 2 ) g ■ fc N 



g-p, 
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ds 



8V27T 2 J ml S-(p- k)'- 



du 



(s,y,P 2 ) [m 2 b - (p + qu) 2 }' 



X- 



2(s-^ 



ds 



) + ( P 2 _ s)u 



8 V2vr 2 [ s - (p - A;) 2 ] 



J^dy(2y-1) 



p2 m c 
2/ 



0,»,p 2 ) mg - (p + gw) 2 



x- 



x(s,y,P 2 
-P 2 



(2q-(p-k)) 2 (l 



2x(s,y,P 2 



u 



g-p' 



(34) 



with the scalar functions $j's and $j's defined by 



(it, f ) = / do; 1 — c<j — i>, u) + 0||(c<j, 1 — oj — t> )) 



duo' 



du' 



1-0)' 



do/' (0 ± (u/', 1 -u/'-u/,u/) + <j>\\(w",l - a;" - a;', a;')) , (35) 



$ 2 (it) 

$ 1 (u,t>) = I dcj ^j_(a;, 1 — — v, v) + 0||(a>, 1 — u — f, v)j 
<£ 2 (it) 



1-oj' 



dw" (0±(u/', 1 - w" - c/, a/) + 0|| (u/ r , 1 - w" - uj\ uj')) . (36) 



The resulting expression (JMj) is much more complicated than those in the case of B — > 
tttt,tcK 17] and B — > J/ibK^La^. because of the mass asymmetry of the two quarks in 
D meson. In contrast to the latter case in which there is no contribution of fa's due to 
cancellation in the corresponding twist-4 parts, the twist-4 piece receive the contributions 
from not only 0j's but also 0j's in the present case. 

To proceed, we should change the above expressions and (f3"4"J) into the desired form 
of double dispersion relation. For the twist-3 contribution F^ k \ we have a dispersion 
expression in (p — k) 2 : 



7 {p-k) 

tw3 



ds 



lm s F t ts k \s, (p + qf,P 2 ,p 2 ) 



IT Jral S — (p — k) 



2 LLLL s J- tw3 



(37) 



where 



lm s F t { ^ k \s, (p + q) 2 ,P 2 ,p 



2 d2 „2\ 



m b f: 



3- 



1QV2 



7T 



dy 



du 



x(s,y,P 2 ) ml — (p + qu) 2 Jx(s,y,P 2 ) V 



^0 37r (l -u,u-v,v) 



X 



s - ^ + tip + q? - P 2 ) (2v - x(s, y, P 2 )) 

y ' 

16^2 



7T 



dy(2y-l) 



du 



dv 



(s,y,P 2 ) mf-(p+ qu) 2 Jx{s,y,P 2 ) V 2 " ^' ^ 



X 



s + — + [(p + q) -p ) (2v-3x(s } y } P 2 )) 

y v ' 



(38) 



10 



Then for Im s F t ^ 3 ^ we make Tailor expansion in variable x(s, y, P 2 ). Up to order 0(x 3 ) the 
result is 



lm s F t ^ k \s, (p + qy,P 
m b f 37T f 1 du 



2 D 2\ 



16\/27r 2 Jo m 2 - (p + qu) 2 



mt. 



s - — + 2v((p + q) 2 -p 2 ) 

y 



u dv 



o ir 



37r (l -u,u-v,v) ((p + q) 2 - p 2 ) 



+ s- 



mz 



(\(p 3n (l -u,u- v,vj] 



y 



+ 



?>7T 



y 



du 



x(s,y,P 2 ) 
x 2 (s,y,P 2 ) 



v=0 



16\/27r 2 Vo m 2 b - (p + qu) 2 

u dv , , 

~2 



1 — U,U — V.v) 



m„ 



-s+ ^ + 2v((p + q) 2 -p 2 ) 

y 



u dv 



o v 



— </> 3 ^(l -u,u-v,v) ((p + g) 2 - p 2 ) 



+ -s + 



+ 



■j(^hn(l-U,U-V,v)^j ^ 

1 - w, u - v , f ) ((p + g) 2 - p 2 ^j 



x(s,y,P 2 ) 



m l\ 9 ( 1 , ,„ 

-S H ^03tt(1 ~U,U-V,V / 

y I ov \v z 



x 2 (s,y,P 2 



0{x 3 



J ii=0 



(39) 



With the substitution u = (m 2 — p 2 )/{s' —p 2 ), the integral in u in the above equation can 
get back to its dispersion form 

ds' F{u{s')) 



(40) 



'o m 2 — {p + qu) 2 Jm\ s' — (p + q) 2 s' — p 2 

At last the desired double dispersion form is achieved. 

The twist-4 contribution in can be treated similarly. The derivation is omitted 

to save some space. We note that compared with the resulting twist-3 contribution, the 
twist-4 part has some additional terms containing denominators of the form 

1 1 



[s — (p — k) 2 ] 



21 2 



or 



[ml — (p + uq)< 



(41) 
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As argued in Q], however, these terms containing higher power of such denominators are 
numerically suppressed and can be neglected safely. Therefore in the ensuing discussion we 
will not take them into account. 

Putting everything together, we obtain the final LCSR result for the soft contribution to 
the matrix element (D°(p)'K (q)\d 1 (0)\B (p + q)): 

—irrib 



(D"(p)^(q)\O 1 (0)\B"(p + q)) i 



SV2^f D f B m 2 B 



1 du 



dy- 



3tt 



— 37r (l -u,u - v,v)\ (2v - x(s,y,m B )) + s 

\ u y 



10 V 

S -H^) ("J_0 37r (i _ u u _ VjV) , ,r(.s. //. /;/7, 

v y J \v 2 J v =o 

i2 



+m b f n 



y 

u dv ~ 
o v 2 



d ( 1 



x 2 {s,y,m 2 B ) 



J D=0 



- u, u - d, v) (3 - y, m B )\ 



( 3 ( \ d 

+ -?0±(1 - u, w - u, u) - -tt'M 1 -u,u-v } v) 



x 2 (s,y,m 2 B ) 



—im b 



v dv 

(l ,/„,K,-*v.^ A 1 ...... 



u=0 



(m 2_ (m 2_ m 2 (1 _ u))/u)/M ,2 



8V27r 2 f D f B m 2 B 

{fs f u dv ( 
~2 [J ^fo«(l-u,u-v,v)^ 



u 



j! A dy{2y-l) 



i\ - 2 



mt — m Df ^ „2- m 



u 



(2v - 3x{s,y,m B )) - s H 

y , 



H ^03tt(1 - w, w - v, v)\ x(s, y, m B ) 



+ 



y J \v 

m 2 b - m 2 D 1 

4 ^03tt(1 -u,u- v,v) 

u v 2 



v=0 



+m b f n 



( m 2 r \ d ( 1 , .„ ,\1 x 2 (s,y,m 2 R ) 

- u, u - d, v) (3 - -ar(s, y, m B )\ 



"(it): 
77 



3 ( -^0±(1 -u,u-v,v] 



-S-<j>±(l -u,u-v,v) 
V ov 



x 2 {s,y,m 2 B ) 



v=0 



,(42) 



where = (m 2 — rri D .)/(s B — m 2 D ). 

Let's proceed to numerical discussion. The D channel parameters are taken as 
m D = 1.87 GeV, m c = 1.3 ± 0.1 GeV, f D = 170 ± 10 MeV ,stf = 6 ± 1 GeV 2 and 
M 2 = 1.5 ± 0.5 GeV 2 . The parameters in B channel are chosen as [23]: tub = 5.28 GeV, 
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m b = 4.7 ± 0.1 GeV, f B = 180 ± 30 GeV, s% = 35 ± 2 GeV 2 and M' 2 = 10 ± 2 GeV 2 . For 
the noirperturbative quantities entering the relevant light-cone distribution amplitudes, we 
use 



22j f 3n = 0.0026 GeV 2 , <5 2 (/i b ) = 0.17 GeV 2 and e{/i b ) = 0.36, whth fi b 



m b /2 ~ 2.4 GeV. With these inputs, the contributions of twist-3 and -4 fall into the following 
ranges respectively: 



i(D°7r |di|5°>^ 3) = (0.024 - 0.053)GeV 3 , (43) 



and 



z^VlOilS )^ = (0.009 - 0.017)GeV 3 , (44) 
the total contribution reads 

i(D -K \6 l \B ) s = (0.033 - 0.070)GeV 3 . (45) 

These sum rule results show a good stability against the variations of both the Borel pa- 
rameters in the given ranges. 



III. DISCUSSION AND CONCLUSION 

Having at hand the LCSR result ()45|) for the matrix element (D ir \Oi\B ) s, we can 
discuss the numerical influence of the power-suppressed soft effect on B° — > D°ti . 

Taking Ci(> 6 ) = -0.257, C 2 (^ 6 ) = 1.117, \V cb \ = 0.043 and \V ud \ = 0.974, the magnitude 
of As given by ((Zj) is estimated at 

\A S {B° — ► D 7r )| = (2.48 - 5.27) x lO^GeV. (46) 

It is in order that we make a numerical comparison between A S (B° — > D°ir°) and the naive 
factorization piece of the decay amplitude Af(B° — > D°ir°) given by ©• With the LCSR 
result F^(m 2 D ) = 0.30 24], we have 

= A S (B° — > D°n°)/A F (B° — D°7T°). 

= 0.54- 1.15. (47) 

This result shows that the resulting soft effect is comparable numerically with the corre- 
sponding factorizable one. A similar ratio was estimated for the B — > J/tpK decay in 



13 



Ref . [bst . with the value 0.30 — 0.70. It seems that power suppressed soft effects are even 
more important in (B° — > D°7r°) than in B —>■ J/tpK, as expected. 

It is also interesting to compare numerically As(B° — > D°ti ) with the factorizable 
contribution of the 2 operator, A { ° 2 \b — ► JJV) = -\C 2 G F V d y* d m\f D F^{m 1 D ). To 
this end, we estimate the ratio 

R 2 = A S (B° — > D\°)/Af 2 \B° — > £>V). (48) 

The result is R 2 = 0.17 — 0.35. Explicitly, our LCSR calculations favor a 2 L > 0, indicating 
that the correction to a 2 , which is relevant to the factorizable and power-suppressed soft 
part, is positive. This forms a striking contrast to the case of Ref. 15[ where R 2 is found 
to be -0.7, a large negative number, so that sign of a 2 L is predicted to be negative, i.e, a 2 
would receive a negative number correction from such non-leading effects. 

Naive factorization does not apply for the color-suppressed B° — ► D°7r° decay and as a 
consequence, the power-suppressed soft exchange correction is expected to be important and 
is worth discussing carefully, in spite of its non-leading character. We discuss such effect in 
the generalized QCD LCSR. The numerical result is in agreement with one's expectations. 
The size of the resulting contribution to the decay amplitude is found to be comparable with 
the corresponding factorizable one, about (50 — 110)% of the latter, and the parameter a 2 
would receive a positive number correction, analogously to the case of B — > J/ipK. These 
observations would be of important phenomenological interests. Of course, at this stage we 
are not able to go a step further to give a complete estimate of the B° — ► D°tt° decay 
amplitude, due to the unknown leading nonfactorizable soft contributions. More theoretical 
or phenomenological efforts in this direction are necessary to better understand the color- 
suppressed charmed decays of B mesons. 
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